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Abstract 

We study the measurement of the atmospheric neutrino oscillation parameters, 623 
and Am23, disappearance channel at three conventional beam facilities, 

the SPL, T 2 K-phase I and NOi^A. These two parameters have been shown to be 
of crucial importance in the measurement of two of the unknowns of the PMNS 
mixing matrix, 813 and the leptonic CP-violating phase 6. In our analyis, the effect 
of the two discrete ambiguities, sign(Am23) and sign(tan 2023), is explicitly taken 
into account. We analyse also the disappearance channel at the Neutrino Factory, 
and combine it with the “golden” iZg —> and “silver” iZg —> appearance channels 
to study its impact on the measurement of 0 i 3 and 5 . Eventually, we present the 
sensitivity of the four facilities to different observables: 0i3, 6, maximal 023, the sign 
of the atmospheric mass difference, Satm, and the 023-octant, Soct- 




1 Introduction 


The results of atmospheric, solar, accelerator and reactor neutrino experiments [1] show 
that flavour mixing occurs not only in the hadronic sector, as it has been known for 
long, but in the leptonic sector as well. The experimental results point to two very 
distinct mass differences^, ~ 8.2 x 10“^ eV^ and ~ 2.5 x 10“^ eV^. 

Only two out of the four parameters of the three-family leptonic mixing matrix Upmns 
[4] are known; 6*12 ~ 34° and 6^23 ~ 45°. The other two parameters, 6*13 and 5, are still 
unknown: for the mixing angle 6^13 direct searches at reactors [5] and three-family global 
analysis of the experimental data [ 6 ] give the upper bound 6*13 < 11.5°, whereas for 
the leptonic CP-violating phase 5 we have no information whatsoever. Two additional 
discrete unknowns are the sign of the atmospheric mass difference and the 023 -octant 
(if 023^45°). 

The full understanding of the leptonic mixing matrix constitntes, together with the 
discrimination of the Dirac/Majorana character of neutrinos and with the measnrement 
of their absolute mass scale, the main neutrino-physics goal for the next decade. In 
the recent past, most of the experimental breakthroughs in nentrino physics have been 
achieved by exploiting the so-called “disappearance channels”: by observing a dehcit in 
the nentrinos that reach the detector with respect to those expected to be emitted from 
the source, a positive and eventnally unambiguous signal of nentrino oscillations has 
been established. The SK detector has gathered indirect evidence of > Vr conversion 
of atmospheric nentrinos. However, no direct observation of r’s is possible at this 
detector. To observe directly fi ^ t conversion, two detectors are under construction 
at the Gran Sasso Laboratory [7,8]. The SNO experiment [9] has shown that a fraction 
of the z/g’s emitted by the Sun core reach the Earth converted into z/^’s and z/t-’s (and 
not into unobservable sterile neutrinos). However, the SNO detector is not able to 
distinguish between r’s and /i’s and thus measure the subleading oscillations z/g —^ 
directly. For this reason, new-generation experiments have been proposed to 
look for the fleeting and intimately related parameters 0 i 3 and 5 through the more 
promising “appearance channels” such as z/g —z/^ or z/^ —z/g (the “golden channel”) 
and z/g —(the “silver channel”). However, strong correlations between 0 i 3 and 5 [10] 
and the presence of parametric degeneracies in the (013,5) parameter space, [11]-[14], 
make the simultaneous measurement of the two variables extremely difficult. A further 
problem arises from our present imprecise knowledge of atmospheric parameters, whose 
uncertainties are far too large to be neglected when looking for snch tiny signals as those 
expected in appearance experiments [15]. Most of proposed solutions to these problems 
suggests the combination of different experiments and facilities, snch as Snper-Beam’s 
(of which T2K [16] is the first approved one), /5-Beam’s [17] or the Neutrino Factory 
[18,19]. 


^ A third mass difference, ~ 1 eV^, suggested by the LSND experiment [2], has 

not being confirmed yet [3] and will not be considered in this paper. 


2 



Clearly, the disappearance channel is the best place to rednce further the uncer¬ 
tainties on atmospheric parameters. It must be reminded, however, that this channel 
is afflicted by a four-fold degeneracy [ 20 ]: the sign of the atmospheric mass difference 
and the 6 ' 23 -octant are poorly measured through this channel, indeed^. It has been 
proposed to face this problem using atmospheric neutrinos at T2K-II, [21,22], or at a 
magnetized iron calorimeter, [23,24], 

The study of the disappearance > z/^ channel at several proposed facilities is the 
main goal of this paper (partially presented in Ref. [25]). How the combination of 
appearance and disappearance channels at different facilities can be used to soften the 
degeneracy problem in the measurement of ^13 and S is also analysed. We first study 
in detail the disappearance channel at three proposed SuperBeam facilities, T2K-I [16] 
(that is under construction at the Jaeri site), NOz/A [26] (that passed first scrutiny 
and will rely upon the existing NuMI facility at FermiLab) and the SPL [27]. We will 
show how energy resolution is crucial to reduce signihcantly atmospheric parameter 
uncertainties and how our present ignorance of ^13 and S affects the results in the 
disappearance channel. 

We then investigate the z/^ —z/^ disappearance channel at the Neutrino Factory. This 
channel, although already considered in the literature (see Refs. [28,29], for example), 
has never gained the center of the stage on its own, being overshadowed by appearance 
channels such as the “golden” z/g —z/^ and the “silver” [30] z/g —^ i^r transitions. This 
channel is able to reduce atmospheric parameter uncertainties to an unprecedented 
level and, using the intertwine with appearance channels, to solve many of the dis¬ 
crete ambiguities for ^13 > 3° — 4°. In this range of di^, combination of appearance 
and disappearance channels is rather effective and 613 and S can also be measured un¬ 
ambiguously. Unfortunately, appearance and disappearance signals are optimized with 
different baselines: degeneracies are solved most effectively with a L = 3000 km base¬ 
line in the appearance channel and with a L = 7000 km baseline in the disappearance 
channel ^ . 

Eventually, we present a comparison of the sensitivity to 6 * 13 , to the sign of the atmo¬ 
spheric mass difference, to the 6 ' 23 -octant and to maximal 623 of the different facilities. 
Their CP-discovery potential [20] is also shown. As expected, the Neutrino Factory 
outperforms the considered SuperBeams in every aspect, with the notable exception of 
the CP-discovery potential that seems to be larger for T2K-II. The Neutrino Factory 
is in this case limited by the “accidental” flow of degeneracies towards |(5| = 0 ° or 180° 
for very small values of 6 'i 3 , see Refs. [33,34]. 

The paper is organized as follows: in Sect. 2 we shortly introduce the four facilities and 
the neutrino-nucleon cross-section; in Sect. 3 we compare the potential of the three 
conventional beams in the measurement of the atmospheric parameters 623 and Am 23 ; 

^ This is why we have to deal with these ambiguities when looking to appearance signals, 
ultimately. 

^ This was known for long (see Refs. [31] and [32] at this regard). 
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in Sect. 4 we study the potential of the Neutrino Factory in the measurement of the 
atmospheric parameters and combine the disappearance channel with “golden” and 
“silver” appearance channels; in Sect. 5 we show the sensitivity to several unknowns of 
the considered facilities combining appearance and disappearance channels; in Sect. 6 
we eventually draw our conclusions. 


2 The experimental setup 

In this section we describe, briefly, the four facilities that we will use in the following 
and we remind the neutrino-nucleon cross-section used throughout the paper. 

2.1 The T2K Beam 

A Super-Beam is a conventional neutrino beam with a proton intensity higher than 
that of existing (or under construction) beams such as K2K [35], NuMI [36] and the 
CNGS [37]. With respect to the /5-Beam [17] and the Neutrino Factory [18], neutrino 
beams of a new design, it has the advantage of a well known technology. On the other 
hand, the flux composition (with as the main component for a tt’*' focusing, plus 
a small but unavoidable admixture of Ve and z/g) limits its sensitivity to Vf. 

oscillations. 

The T2K-I facility [16] has been approved and it is under construction. It consists of 
a conventional beam with 0.75 MWatt power 50 GeV protons produced at the J-Parc 
site in Tokai aiming at the Super-Kamiokande water Cerenkov detector with a baseline 
L = 295 km. The characteristic feature of the T2K beam is (apart from its power, 
that makes this beam a first-generation Super-Beam) that the detector is not on the 
longitudinal axis of the beam. This feature has a threefold advantage: first, it allows 
that neutrinos produced through protons whose energy is fixed by the requirement 
that they may be used for different purposes have a z/e oscillation probability 

peaked at L = 295 km, where the detector is located; second, the off-axis beam is 
narrower than the on-axis one, thus improving the matching of the LjE ratio with the 
hrst peak of the oscillation; third, it allows a reduction of the beam-driven background. 
Notice that the off-axis angle has not yet been chosen, to arrange the LjE ratio to 
slight variations in the measured value of The technical design of the tunnel 

is such that the off-axis angle can be modihed ranging between 2° and 3°. The T2K 
fluxes for a 2° off-axis angle as they are expected at the Kamiokande site are shown in 
Fig. l(left). Following the Letter of Intent [16], we have considered four 200 MeV bins 
with Ejnin = 400 MeV. The muon and electron identihcation efficiencies per bin are 
presented in Fig. 21. The average neutrino energy is < E^, >= 0.75 GeV with 5 years 
of TT’*' running, only. 

The dominant background is made of pions from neutral currents. The detector is ex¬ 
tremely well understood, being the SuperKamiokande water Cerenkov. As for the SPL, 
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Fig. 1. Left: T2K-1 fluxes at the Kamioka location (295 km baseline), [16]; Right: NOuA 
fluxes at the far location (812 km baseline), [38]. 

the dominant sources of systematic error (apart for the neutrino-nucleon cross-sections) 
are the dehnition of the hducial volume in the detector and the overall normalization 
of the flux. We have introduced a “pessimistic” 5 % global systematic error. 

A second phase has been proposed for the T2K experiment: if a positive signal in 
conversion is observed, the beam power is planned to be increased up to 4 
MWatt and the detector mass up to 1 Mton (500 Kt hducial mass). A 5 years vr"'' and 
5 years vr" running are envisaged. Such scenario will be called T2K-II in Sect. 5. 

2.2 The NOvA Beam 

The NOz/A facility Letter of Intent [26] has passed the hrst scrutiny. It consists of a 
conventional beam with 0.7 MWatt power 120 GeV protons produced at the NuMI in 
FermiLab aiming at a site located at L = 812 km from the source. It is still not decided 
the typology of the detector that will be located there, with two options considered 
so far: a low-Z calorimeter; or, an Icarus-like liquid argon bubble chamber. As for 
T2K, the NOz/A detector will be placed at a non-zero off-axis angle with respect to 
the longitudinal beam direction, to reduce backgrounds and to have a narrow beam to 
match the L/E oscillation peak ratio. A reduced version of the detector will be placed 
at the FermiLab site to extrapolate the expected rates and their energy spectrum in 
the far detector in the absence of oscillation. The NOuA fluxes for a 0.72° off-axis angle 
as they are expected at the far site are shown in Fig. 1 (right). 

In the following, we consider a totally active 30 kton detector made of liquid scintillator 
and PVC [40] with 6.5 x lO^*’ pot/y (no proton driver is considered, [39]). The average 
neutrino energy is < E^, >= 2.22 GeV, with 5 years of 7 i~^ run only [26]. Events are 
grouped in three ~ 650 MeV bins with Emin = 1000 MeV, with a constant muon 
identihcation efficiency = 0.9 [41] and a constant electron identihcation efficiency 
Ce = 0.24 [40], 
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The dominant background for the disappearance channel are Vx neutral currents 
with a rejection factor at the level of 5 x 10“^ (for the appearance channel the back¬ 
ground is typically from electron neutrinos in the beam and from neutral current events 
faking electron neutrinos, with a rejection factor approximately of 2 x 10“^). 

The main source of systematic errors are the incomplete knowledge of the neutrino 
flux and the extrapolation of the flux and of the backgrounds from the near to the far 
detector. The near detector will be used to improve the present knowledge of low energy 
neutrino cross-sections and to complement the results of the MINERz/A experiment 
[42] (operating before the bulk of NOz/A data is collected). We have introduced a 
“pessimistic” 5% global systematic error. 


2.3 The SPL 

The SPL is a proposal for a high intensity conventional beam with a 4 MWatt 2.2 GeV 
proton driver located at CERN, Ref. [27]. The neutrino fluxes have first been computed 
in a full simulation of the beamline in Ref. [43], assuming a decay tunnel length of 60 m. 
The corresponding fluxes are shown in Fig. 2(left). However, the original beam design 
was conceived as the first stage of a Neutrino Factory, and it was not optimized as 
a facility to look for z/g oscillations on its own. Such an optimization has been 

presented in Ref. [44], and the modihed fluxes are shown in Fig. 2(right). We will 
consider both the “old” and “new” neutrino fluxes to study the potential of the z/^ 
disappearance channel at this facility to measure the atmospheric parameters 023 and 

For both fluxes, a UNO-like [45] 1 Mton water Cerenkov detector (with a 440 kton fidu¬ 
cial volume) located inside the Frejus tunnel with a baseline L = 130 km is considered. 
The average (anti)neutrino energy for both the “old” and the “new” fluxes is 

Ey G [250 — 300] MeV, such as to match the first oscillation peak. 2 years of running 
and 8 years of 7r“ running are considered. The detector characteristics when exposed 
to the SPL beam were studied in detail in Ref. [46], where an almost constant Cg = 70 
% efficiency to identify electrons was found. Reasonably enough, due to the low energy 
of the neutrino flux, no binning was considered (being the neutrino energy comparable 
with the Fermi motion). Following Ref. [47], however, we try to consider a very coarse 
binning grouping events in two bins with energies Ei G [0, 250] MeV and E2 G [250, 600] 
MeV. We have considered a constant bin independent muon identihcation efficiency, 
= 0.7. The electron identihcation efficiency is Cg- = 0.707, eg+ = 0.671 [48]. 
Notice that in the appearance channel no binning is considered, since in the energy 
range Ei practically no event is observed. 

The background is made dominantly of pions from neutral currents and it is negligible 
in the range of energy E 2 . The dominant sources of systematic error (apart for the 
neutrino-nucleon cross-sections) are the dehnition of the hducial volume in the detector 
and the overall normalization of the hux. We have introduced a “pessimistic” 5 % global 
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Fig. 2. Left: Standard SPL Super-Beam fluxes at the Frejus location (130 km baseline), [43]; 
Right: Optimized SPL Super-Beam fluxes at the Frejus location (same baseline) [44]- 

systematic error. 


2-4 The Neutrino Factory 

The Neutrino Factory proposal considered here is a 50 GeV muon storage ring fuelled 
through a 4 MWatt SPL-like Super-Beam [49], with 2 x 10^° muons decaying in each 
straight section of the storage ring per year. Five years of data taking for each muon 
polarity are envisaged. Two detectors of different technology are considered: a 40 kton 
Magnetized Iron Detector (MID) located at L = 3000 km or L = 7000 km; and a 
10 kton Emulsion Cloud Chamber (ECC) at L = 732 km. This proposal corresponds 
to the design of a possible CERN-based Neutrino Factory Complex, with detectors 
located at the Gran Sasso Laboratory as the shortest baseline and at a second site 
to be dehned as the medium and long baseline. These detectors have been especially 
optimized to look for a particular signal: the MID for the “golden” channel z/g — 
and the ECC for the “silver” channel zzg ^ 1 /^-. The corresponding neutrino fluxes are 
shown in Fig. 3(left). 

The detectors background, muon identihcation efficiency and systematics for the golden 
and silver channels at this specific facility have been studied in Refs. [50] (the Magne¬ 
tized Iron Detector) and [51] (the Emulsion Cloud Chamber). It has been noticed that 
the MID simulation must be updated due to the growing evidence for the relevance of 
low energy bins that were sacrified at first to reduce backgrounds in Ref. [10]. To take 
advantage of the low energy bins, in this paper we will consider for the disappearance 
channel 4 GeV bins when the MID is located at the L = 3000 km baseline and 5 GeV 
bins when the MID is located at the L = 7000 km baseline, instead of the 10 GeV 
bins considered in the appearance channels. The disappearance signal is much differ¬ 
ent from the golden and silver appearance signals. In particular, not such strong cuts 
as for the appearance channels signals must be applied to reduce the background. A 
reasonable constant muon identification efficiency = 0.7 has been applied when the 
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Fig. 3. Left: 50 GeV Neutrino Factory fluxes at 3000 km [10]; Right: the vN and vN 
cross-sections on water [53]. 

disappearance channel is considered throughout the paper. A global 2 % systematic 
error has been included. 


2.5 The neutrino cross-seetion 


An important source of systematic error is our present poor knowledge of the vN and 
vN cross-sections for energies below 1 GeV [52]: either there are very few data (the 
case of neutrinos) or there are no data at all (the case of antineutrinos). On top of 
that, the few available data have generally not been taken on the target used in the 
experiments (either water, iron, lead or plastics), and the extrapolation from different 
nuclei is complicated by nuclear effects that at the considered energies can play an 
important role. 

We use for each detector a different cross-section following Ref. [53]. For example, we 
show in Fig. 3(right) the cross-sections on water used for the water Cerenkov detectors 
throughout the paper. Notice the difference between the V(,N and cross-sections: 
the former, being an interaction between the Ve and a neutron inside the oxygen nucleus, 
is affected by nuclear effects and thus shows a threshold energy. The latter is mainly a 
Vf. interaction with the protons of the two hydrogens, approximately free. This effect, 
although less pronounced, is visible also for and for 0 ^. This feature is quite relevant 
for neutrino/antineutrino of hundreds of MeV energy, region where different cross- 
sections can easily differ by a factor 2. Be aware that there are other nuclear effects 
(see Ref. [54] and references therein) not included yet in any of the available calculations 
that could play an important effect at the cross-section threshold energy. 

In the rest of the paper, we will make use of cross-sections on iron, lead and scintillator 
for the MID, the ECC and the NOVA detectors, respectively. 
































3 The measurement of 623 and at SPL, T2K and NOz/A. 


3.1 Correlations and degeneracies in disappearance 

A conventional beam with neutrinos of moderate energy can perform an independent 
measurement of the atmospheric parameters 623 and Am 23 disappearance 

channel; it is expected that this kind of facilities will reduce the error on the atmospheric 
mass difference to less than 10 % in a few years for > 2.2 x 10“^ eV^ [55]. The 

expected error on the atmospheric angle depends on the value of 623 itself, the smallest 
error for large but non maximal mixing (as it has been shown in Ref. [56]). It is 
interesting to study in detail the parameter correlations and degeneracies that affect 
this measurement and that can induce large errors. We remind the vacuum oscillation 
probability expanded to the second order in the small parameters 613 and {A 12 L/E) 
|57]: 



where J = cos ^13 sin 2^12 sin 2^13 sin 2^23 and A 23 = Am23/‘2E, A 12 = Ami 2 /‘^E- The 
hrst term in the first parenthesis is the dominant one and is symmetric under 623 
7 r /2 — 6 * 23 . This is indeed the source of our present ignorance on the 6 ' 23 -octant [13], 
parametrized by Soct = sign(tan2023)- This symmetry is lifted by the other terms, that 
introduce a mild CP-conserving ^-dependence also, albeit through subleading effects 
very difficult to isolate. 

Considering that the sign of the atmospheric mass difference, Satm = sign(Am 23 ), is 
unknown at present, for any experimental input { 623 , we must solve two systems 

of equations: 

^i(^23, Aml^^] Satm) = N^f,{023, \ Aml^l] S atm) , (2) 

^i(^23, Amlt^] Satm) = N^f,{023, lAm^g]; -Satm) , (3) 


where Satm is the physical mass hierarchy. 

For non maximal 623 four different solutions are found: for |Am 23 | ^^Itm we get 
two solutions from eq. (2), i.e. the input value 623 = O 23 and 623 ~ 7 r /2 — 623 , being the 
second solution not exactly at 623 = 71/2 — 623 due to the small 6 ' 23 -octant asymmetry; 
and two more solutions from eq. (3) at a different value of lAm^gl [20]. In eq. (1) 
we can see that changing sign to Am^s the second term becomes positive: a change 


9 



3 



2 I-.-.-^- 

0.3 0.4 0.5 0.6 0.7 

sin^ 023 , 


3 



2 I-^^-.- 

0.3 0.4 0.5 0.6 0.7 

sin^ 023 


Fig. 4. The sign degeneracy at T2K-I; left: 023 = 45°; right: 023 = 41.5°. 
that must be compensated with an increase in |Am 23 | to give Satm) = 

The result of a fit to the disappearance channel data at the T2K-phase I experiment 
is depicted in Fig. 4 for three different values of the atmospheric mass difference: 
Am 23 = (2-2, 2.5, 2.8) x 10“^ eV^. Fixed values of the solar parameters have been used, 
Am 32 = 8.2 X 10“® eV^; 0i2 = 33°. For maximal atmospheric mixing 023 = 45°, Fig. 4 
(left), two solutions are found at 90 % CL when both choices of Satm are considered. On 
the other hand, using a non maximal input atmospheric angle 023 = 41.5° (sin^ 023 = 
0.44) [6] four degenerate solutions are found. Fig. 4(right). In general, we must therefore 
speak of a two-fold or four-fold degeneracy in the disappearance channel, as it was 
pointed out in Ref. [20]. 

Notice how the disappearance sign clones appears at a value of |Am 23 | higher than 
the input value, as it was expected from eq. (1). The shift in the vertical axis is a 
function of 0 i 3 and 6 , that in this particular case have been kept hxed to 0 i 3 = 0° = 5. 
The degeneracy can be softened or solved by using detectors with large baselines that 
can exploit matter effects, as it will be shown in Sect. 4. Notice, eventually, that the 
uncertainty on the atmospheric parameters can be enhanced once we take into account 
that 013 and 6 are unknown [56]. This will be studied in Sect. 3.4. 

3.2 A matter of conventions 

It is useful to open here a short parenthesis to address a problem that arised re¬ 
cently concerning the “physical” meaning of the variable used to £t the “atmospheric” 
mass difference, Notice, first of all, that the experimentally measured solar 

mass difference Am^^; can be unambiguously identified with the three-family param¬ 
eter Am 32 = ~ This is not true for the experimentally measured atmospheric 

mass difference Since the subleading solar effects are, at present, barely seen in 
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atmospheric neutrino experiments, we can define in different ways the three-family pa¬ 
rameter to be used in the fits: for example, using = m^—rn^ (the choice adopted 

throughout this paper), or Am ^3 = — m\, or even Aw? = (Am^s + Am?3)/2, [58], 

we will get in general the same results using present data. When future experiments 
aiming to the measurement of the atmospheric mass difference at the level of 10“"^ eV^ 
will be running, however, different choices of the fitting parameter will give different 
results. 

This can be observed in Fig. 5, where the three choices introduced above are compared. 
We plot in the three panels the 90 % CL contours resulting from a fit to the experimental 
data corresponding to the input value, Am^^^ = 2.5 x 10“^, in normal hierarchy, 
but fitted in turn in Am"^^ (left panel), Am^g (middle panel) and Am^ (right panel). 
As it can be seen in the figure, the contour corresponding to the normal hierarchy, 
Satm = Satmi IS always locatcd around the input value (in each plot the input value 
corresponds to a different fitting variable, though). On the other hand, the contour 
obtained for the inverted hierarchy is located above, below and on top of the input 
value, respectively, depending on the fitting variable. This is a clear consequence of 
the fact that the difference between each of the possible choices is 0 {Aml 2 ) and it is 
reflected in a different form of eq. (1). 

Up to here, it is perfectly clear what happens whenever we use a certain three-family 
variable to fit the results that at present are given under the label A philo¬ 

sophical discussion, however, arised around the “physical” meaning of the different 
choices reported above. The idea is that the “physically meaningful” quantities to be 
measured in oscillation experiments are the “oscillation frequencies”. For three-family 
mixing, three “frequencies” can be defined, the shortest being the solar oscillation fre¬ 
quency (as we said, unambiguously related to the mass difference Am‘l 2 )- In normal 
hierarchy, the middle frequency is related to Am^g and the longest one to Am\^. In 
inverted hierarchy these two frequencies are interchanged and the middle frequency will 
be related to Am^ and not to Am| 2 - For this reason, it has been suggested that plots 
in different hierarchies should be presented using different variables where maintaining 
the ordering of the oscillation frequencies. For example, if we choose to identify Am^^^ 
with Am 23 normal hierarchy (i.e., with the middle frequency), we should identify 
it with Am ‘^1 in inverted hierarchy (i.e., again with the middle frequency). As a con¬ 
sequence, we could not plot contours for the two hierarchies on the same figure (the 
vertical axis corresponds to different variables depending on the choice of Satm)- This is 
a drawback of giving to the oscillation frequencies a deeper physical meaning than to 
the mass differences. A new variable to be identified with Am\^^ has been introduced 
in Ref. [58]: Am^ = (Am^g -|-Am^g) /2. When changing hierarchy, this variable just 
flip its sign (since Am^g —> — Amfg and viceversa). It is thus possible to present both 
hierarchies in the same figure with |Am^| in the vertical axis. The oscillation frequency 
related to this variable is neither the longest nor the middle one, but it maintains its 
role of next-to-longest or a-bit-longer-than-the-middle-one in both hierarchies ^ . 

^ It should be stressed that even fitting in inverted hierarchy using Am^ as the fitting vari- 
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Fig. 5. Different choices of the three-family “atmospheric” mass difference; left: Am 23 ; ™ddle: 
Aml^; right: Am?, [58]. 

It seems to us that the “physical meaning” of a frequency is not deeper than that of a 
mass difference, and that it is perfectly acceptable to use any variable to £t the exper¬ 
imental data when considering a full three-family analysis. What is really important 
is to be consistent with the adopted choice, in particular when adding appearance and 
disappearance data, something that we will do in Sect. 3.4. In the rest of the paper we 
adopt Am 23 as fitting variable. 

3.3 The importance of energy resolution 

It is extremely important that an experiment whose goal is to improve significantly the 
present uncertainties on the atmospheric parameters may be able to use energy depen¬ 
dence. A counting experiment is certainly limited, as it has been shown in Refs. [59,20]. 

In Fig. 6 we present a comparison of the disappearance channel at T2K-phase I (left 
panels) and NOz/A (right panels). Both maximal mixing, 023 = 45° (top panels), and 
non maximal mixing, 023 = 41.5° (bottom panels), are used as input values. In each 
plot we present the 90% CL contours in the (023,Am23) plane for different energy 
bins (only the bins corresponding to neutrinos with energy just below and above the 
peak energy are reported for T2K, all bins for NOz/A). Again, solar parameters are 
kept fixed to their present best fit values, Am ^2 = 8-2 x 10“® eV^, 0i2 = 33°, and the 
unknowns are hxed to 0i3 = 5 = 0°. Both for T2K-phase I and NOz/A 5 years of 
are considered. Notice that the contours corresponding to neutrinos with an energy 
below and above the oscillation peak have a different shape. As a consequence, the 
combination of different bins signihcantly increases the 023 resolution of the experiment 
with respect to a counting experiment. It can be seen that T2K-I is able to measure 
Am 23 with a precision of less than 10“^ eV^ for Am^g = 2.5 x 10“^ eV^. However, it 
is not able to exclude maximal mixing at 90 % CL for 023 = 41.5°. The same happens 
for NOz/A, see Fig. 6(right). 

In Fig. 7 we present a comparison of the disappearance channel at the SPL using 
the neutrino fluxes of Ref. [43], Fig. 2(left panels), or the neutrino fluxes of Ref. [44], 

able, the allowed region is shifted with respect to that corresponding to the normal hierarchy. 
This can be observed in Fig. 5(right), and it implies that \Am?\i^H A \Am?\iH. 
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Fig. 6. Binning at T2K-I (left) and NOvA (right); top: O 23 = 45°; bottom: 623 = 41.5°. 

Fig. 2(right panels), again for both maximal mixing, 623 = 45° (top panels), and non 
maximal mixing, 623 = 41.5° (bottom panels). The 90 % CL contonrs corresponding to 
both SPL energy bins are drawn. In this case, both nentrinos and antinentrinos flnxes 
are prodnced by tt"'' and 7 i~ decays with 2 and 8 years of data taking, respectively. 
Notice that, however, being the nentrino and antinentrino average energies extremely 
similar in this setnp, the nentrino (antinentrino) contonrs are almost snperimposed. As 
a conseqnence, experimental information from the nentrino and the antinentrino flnxes 
is not complementary and we observe jnst an increase in the statistics. For this reason, 
althongh in this setnp a 1 Mton detector is considered, the resolntion in 623 is not 
astonishing. A big improvement with respect to the resnlts presented for this facility in 
Ref. [20,15] is represented by the spectral information: in the new analysis two energy 
bins are considered for both setnps. Thanks to this, maximal mixing can be exclnded 
at 90 % CL for 623 = 41.5°. As a hnal comment, notice that the new nentrino flnxes 
of Ref. [44] do not improve the 623 resolntion with respect to old flnxes from Ref. [43]. 
This is becanse the new flnxes have been optimized to look for the —>■ appearance 
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Fig. 7. Binning at the standard (left) and optimized SPL (right); top: 623 = 45°; bottom: 
023 = 41.5°. 

signal and not to the —>■ disappearance one. In particular, the average neutrino 
and antineutrino energies are identical (see Fig. 2) and thus the small complementarity 
of the two fluxes reduces. In the rest of the paper we will consider the SPL with old 
fluxes, only. 

It is remarkable how a relatively “small” experiment such as T2K-I, with only 5 years 
of data taking in one polarity and with a 22.5 Kton water Cerenkov detector, has a 
resolution in and 023 not much worse than the SPL facility, with a much larger 

1 Mton water Cerenkov and 10 years of data taking in both pion polarities (compare 
the left panels of Figs. 6 and 7). One of the reasons is that the SPL has both neutrino 
and antineutrino beams with an average energy corresponding to the oscillation peak 
for the L = 130 baseline. As a consequence, information coming from the two beams 
just add statistically but it is not complementary (in the absence of matter effects). It 
would be better to run the SPL with neutrinos again but at a different energy (not at 
the oscillation peak) after the hrst two years with tt’*'. 
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Fig. 8 . The impact 0/^13 at T2K-I; left: O 23 = 45°; right: O 23 = 41.5°. 

3.f The impact 0 / 6*13 ^ 

Up to this moment we have kept the unknown parameters ^13 and 5 as external hxed 
quantities, ^13 = 5 = 0 °, following what we have done for the solar parameters Am ^2 
and 6 * 12 . However, the two sets of parameters should be treated differently. Indeed, we 
do have a good measure of solar parameters and it has been shown in the literature 
that the impact of solar parameter uncertainties in the measurement of atmospheric 
parameters is negligible [15]. The main effect is the shift in the atmospheric mass 
difference htting variable, as it has been discussed in Sect. 3.2. This is not the case 
for the unknown parameters (6*13, d). Being both parameters unknown, we must £t the 
atmospheric parameters introducing them as free variables to be reconstructed at the 
same time with (^ 23 , ^ 77 x 23 ). ^ recent comprehensive three-families analysis of present 
solar, atmoshperic and reactor data can be found, for example, in Ref. [ 6 ]. 

In Fig. 8 we show the effect of a varying ^13 (but hxed S) on the previous hts at 
T2K-I for 6*23 = 45° (left panel) and 023 = 41.5° (right panel) and three values of the 
atmospheric mass difference, Am 23 ~ (2.2, 2.5, 2.8) x 10“^ eV^. Both choices of Satm 
are shown on the same plot. The input values of the unknowns are: ^13 = 0°, 5 = 0°. 
In the £t, 6*13 is free to vary in the range 6*13 G [0°, 10°]. Notice that the input values 
can be htted with increasing values of 6*13 if 6*23 is also increased, resulting in a shift of 
the 90 % CL contours to the right. This is a consequence of the 6*23- asymmetric second 
term in the hrst parenthesis of eq. ( 1 ). 

Such a naive treatment of the experimental data is, however, not correct. The con¬ 
sidered Super-Beam facilities have indeed been proposed to look for the appearance 
t'/i channel, the oscillation probability in vacuum expanded to the second order 

in the small parameters 6*13 and {^i 2 L/E) [10,57] for which is: 

Ue) = §23 sin^( 26 * 13 ) sin^ (“f“) 
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+ sin(26'i2) sin( 26 ' 23 ) sin(26'i3) cos sin (~^) 

+ C23 sin2(20i2) sin^ (~f“) ’ 

where ± refers to neutrinos and antineutrinos, respectively. It is clear that to take 
properly into account the effect of ( 6 * 13 , 5) we should combine (whenever possible) ap¬ 
pearance and disappearance signals at a given facility. 

The results of a four-parameters fit in ( 6 * 23 , 6 * 13 , 5) projected onto the ( 6 * 23 , Am^a) 

plane obtained by combining the disappearance and appearance signals at the T2K-I, 
NOz/A and SPL facilities are presented in Figs. 9-11, respectively. In all hgures two 
choices of 623 and ^13 are shown, 623 = 45°, 41.5° (left and right panels) and ^13 = 0°, 8 ° 
(top and bottom panels), whereas 5 = 0°. Solid lines represent the result of a £t with 
variable 6*13 and 6 . As a reference, we also present the results of a £t ® with 613 = 613 
and 6 = 6 (dotted lines). It can be observed that in most of the cases the combination 
of disappearance with appearance signals reduces the spread in 623 that was observed 
in Fig. 8 . Only for ^13 large we can still see some effect. 

Eventually, the effect of a non-vanishing 6 can be easily understood: the CP-conserving 
5-dependent term in the second row of eq. ( 1 ) is positive for |5| G [0°,90°], negative 
for |(5| G [90°, 180°] and vanishing for |5| = 90°. When non-vanishing, this term induces 
a shift in the htted value of Am^s that is maximal (minimal) for 5 = 0 ° (5 = 180° ). 
Being, however, a 6 ' 13 -suppressed perturbation on the subdominant 0 {Aml 2 ) term, its 
effect is rather small and it will not be shown in a separate figure. 

As a hnal comment, we notice that the three considered Super-Beams are indeed ca¬ 
pable of strongly reducing the present uncertainties on the atmospheric parameters 
623 and Am 23 , particular for large values of ^13 when combination of appearance 
and disappearance signals is most effective. For 0i3 = 8 °, for example, maximal atmo¬ 
spheric mixing can be excluded at 90 % CL at all three facilities for 623 = 41.5° (only 
the SPL, with its gigantic 1 Mton water Cerenkov, is able to exclude maximal mixing 
for Q 23 = 41.5° for a vanishing ^ 13 ). Notice, however, that while the disappearance oc¬ 
tant degeneracy appears to be solved when combining disappearance and appearance 
signals for ^13 = 8 ° when ^13 is treated as a fixed parameter, is indeed restored if 613 is 
left freely varying in the presently allowed range. 


4 Appearance and disappearance channels at the Neutrino Factory. 

The disappearance signal has not yet been deeply investigated in the framework of the 
Neutrino Factory (see, however. Ref. [29]). The top-charted signals have been up to now 
the “golden” z/g —z/^ [10] and the “silver” z/g — r'r [30] channels. However, the Neutrino 
Factory produces a very intense z/^ beam that can be used to measure z/^ —z/^ and the 

^ Only the results for Satm = Satm are shown. 
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Fig. 9. Appearance and disappearance at T2K-I. Left: 623 = 45°; right: 623 = 41.5°; top: 
9 13 = 0°; bottom: 613 = 8°. 

leading > Vj. appearance channel as well. This last transition is probably the best 
signal to improve signihcantly the atmospheric parameter nncertainties. However, to 
stndy it in detail we need a detector able to identify efficiently an enormons amonnt of 
r’s (ronghly 1000 and 400 r“ and are expected, respectively), snch as it is expected 
in a Nentrino Factory beam. This detector can be in principle the same used to look 
for the silver channel, but the needed scaling of the scanning power must be carefully 
studied and it will not be discussed here. We will devote this Section to a detailed 
study of the combination of golden, silver and disappearance channels at the Neutrino 
Factory. 

Two possible baselines are studied, L = 3000 km and L = 7000 km, i.e. the optimal 
distance to look for a CP-violating signal [10] and the optimal distance to exploit 
matter effects through the disappearance channel, respectively. 

In Figs. 12 and 13 we hrst present the 90 % contours of the on-peak and above-peak 
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Fig. 10. Appearance and disappearance at NOuA. Left: 623 = 45°; right: 623 = 41.5°; top: 
013 = 0 °; bottom: 0 i 3 = 8 °. 

disappearance channel energy bins for two valnes of 023, 023 = 41.5°, 45° (left and 
right panels, respectively), and two valnes of 0 i 3 , 0 i 3 = 0 °, 8 ° (top and bottom panels, 
respectively). The medinm baseline resnlts are shown in Fig. 12, the long baseline 
resnlts in Fig. 13. Both nentrino and antinentrino bins are presented. Again, solar 
parameters are kept fixed to their present best fit valnes, Am ^2 = 8-2 x 10 “^ eV^, 
012 = 33°. The inpnt valne for h is h = 0°. 

Notice, first of all, that the resolntion in 023 is extremely good and that maximal mixing 
can be easily exclnded for 623 = 41.5°. In particnlar, at the L = 7000 baseline a 4 % 
error on 623 for 623 = 41.5° is fonnd. This was expected, being the statistics mnch higher 
than at the Snper-Beams experiments stndied in Sect. 3 (see Tabs. 1-5). Moreover, a 
new featnre arises for 0 i 3 = 8° at both baselines. When a rather large non-vanishing 
013 is switched on, matter effects become extremely important and introdnce a strong 
023-asymmetry in eq. (1). The asymmetry can be clearly seen in the bottom panels of 
Figs. 12 and 13, and it is crncial in solving the disappearance oetant degeneraey, see 
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sin^ @23 sin^ 623 

Fig. 11. Appearance and disappearance at the standard SPL. Left: 623 = 45°; right: 
023 = 41.5°; top: 613 = 0°; bottom: @13 = 8 °. 

Sect. 3.1. 

We point out that the on-peak bin shows a circular shape centered around the input 
value, distinct from the upward(downward)-curved shapes of contours corresponding 
to above(below)-peak bins of Figs. 6 and 7. This is because the Neutrino Factory flux 
is not centered around the peak energy for the chosen baselines: the peak energy for 
a L = 3000 km baseline would be < >~ 6 GeV and for a L = 7000 km baseline 

< >~ 14 GeV. As a consequence, we have no energy bins below the peak energy, 

but only on-peak or above-peak bins. This is a major flaw of the present Neutrino 
Factory design, that could perhaps be solved with an improved detector capable to 
take advantage of low energy bins. 

The solving of the octant-degeneracy in the disappearance channel (for large ^ 13 ) is of 
great importance. It is useful to recall that, when looking for a GP-violating signal in 
appearance channels as z/g —>■ z/^, z/g —^ Vj. and ^ Vei we must deal in general with an 
eightfold-degeneracy [14] that originates from three sources: the discrete ambiguities 


19 




























sin^ 023 
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Fig. 12. Binning at the L = 3000 km Neutrino Factory. Left: 023 = 45°; right: 023 = 41.5°; 
top: 013 = 0°; bottom: 0 i 3 = 8°. 

parametrized by Satm [12] and Soct [13] and the intrinsic ambignity due to the trigono¬ 
metric nature of Pe^i, Per and P^e in the two unknowns 0i3 and S [11] (see eq. (4) as an 
example). As a consequence, for generic values of the input pair (0i3,h), we get eight 
different solutions in the (0i3,5) plane. Different ways to solve the eightfold ambiguity 
have been proposed in the literature, such as to combine experiments with different 
LjE [11], or the golden channel at the Neutrino Factory with —>■ z/g at a Super- 
Beam [33], or the golden and silver channels at the Neutrino Factory [30,51], or even to 
combine different SuperBeams [60,61] or a SuperBeam and a /d-Beam [62,59] (see also 
Refs. [63]-[67]). In Ref. [68] it was shown that the combination of the Neutrino Factory 
with a SPL-like SuperBeam was able to solve the eightfold-degeneracy and to give a 
single allowed region in the ( 0 i 3 , 5) plane at the price of three specialized detectors, 
a 40 kton Magnetized Iron Calorimeter and a 4 kton Emulsion Cloud Chamber (to 
look for z/g —>■ Vfj, and z/g —^ Vj. oscillations, respectively) and a gigantic 1 Mton Water 
Cerenkov (to look for z/^ ^ z/g). The results of Figs. 12 and 13 suggest, on the other 
hand, that some reduction of the eightfold ambiguity can be achieved at the Neutrino 
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Fig. 13. Binning at the L = 7000 km Neutrino Factory. Left: 023 = 45°; right: 023 = 41.5°; 
top: 013 = 0°; bottom: 0 i 3 = 8°. 

Factory combining both golden and silver appearance signal with the disappearance 
signal, with no need for a huge water Cerenkov detector to be located somewhere else 
(albeit, using a 10 kton ECC). 

The results of a four-parameters £t in (023, Am 23 , 0i3, 5) projected onto the (023, Am^g) 
plane® obtained by combining the disappearance and appearance signals at the Neu¬ 
trino Factory are eventually presented in Fig. 14 for the medium baseline L = 3000 
km and in Fig. 15 for the long baseline, L = 7000 km. Again, two values of 023, 
023 = 41.5°, 45° (left and right panels, respectively) and two values of 0i3, 0i3 = 0°,8° 
(top and bottom panels, respectively) are shown. 

The input value for 5 is 5 = 0°. In each figure, dotted lines refer to the fit with 0i3 = 0i3, 
5 = 5 and Satm = Satm- The solid lines, on the other hand, represent the result of a 
£t with variable 0i3,5 and Satm- Notice that, as a general result, the uncertainties 


® The two-parameters 90 % CL contours are shown. 
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Fig. 14. Appearance and disappearance in the ( 623 , P^CLn-e at the L = 3000 km Neutrino 

Factory. Left: 623 = 45°; right: 623 = 41.5°; top: 6*13 = 0°; bottom: ^13 = 8°. 

on the atmospheric parameters are better measured at the long baseline than at the 
medium one. In particular, for ^13 = 0 °, maximal mixing cannot be excluded for a low 
atmospheric mass difference, = 2.2 x 10 “^ eV^ at the medium baseline: for such 

a low value of Am 23 , the energy bins are above-peak and thus show an upward- 
curved shape (see Fig. 12), reducing their complementarity. For 6^13 = 8 ° the octant 
ambiguity is always solved by the strong matter effects at both baselines. Combining 
the disappearance channel with the long-studied golden and silver appearance channels 
at the Neutrino Factory solves the octant-ambiguity for ^13 > 3°. 

In Fig. 16 we eventually show the results of a four-parameters £t in (^ 23 , Am 23 , dl3, (^) 
projected onto the “classic” ( 6 * 13 , 5) plane [ 10 ] obtained by combining the disappearance 
and appearance signals at the Neutrino Factory, both for the L = 3000 km and L = 
7000 km baselines (left and right panels, respectively). Two values of 5 are shown, 
5 = 45° (top panels) and 5 = —90° (bottom panels), and three values of ^ 13 , ^13 = 2°, 5° 
and 8° (corresponding to sin 6*13 = 0.04,0.09,0.14). Solid lines correspond to 90% CL 
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Fig. 15. Appearance and disappearance in the ( 623 , at the L = 7000 km Neutrino 

Factory. Left: 623 = 45°; right: 623 = 41.5°; top: 6*13 = 0°; bottom: ^13 = 8 °. 


contours with {Satm.Soct) = {Satm.Soct), daslicd lines to {Satm.Soct) 7 ^ iSatm.Soct)- As it 
has been known for long, the Neutrino Factory experiments have such a long baseline 
that the sign of the atmospheric mass difference is measured down to 6*13 = 2°. The 
novelty here, is in that at 90 % CL not only the sign clones but the octant and intrinsic 
clones are also solved down to 6*13 = 2° for most values of 5. This stands true for both 
baselines. However, it is manifest that for L = 7000 km the sensitivity to 5 is lost (this 
was called the “magic baseline” [32] precisely because around this distance the effect 
of 5 in the —>■ oscillation probability vanishes, see also Refs. [31] and [14]). A 
very good measurement of 6^13 can be performed at this baseline, but certainly not of 
the leptonic CP-violating phase. The case of <5 = —90° is particularly bad, due to low 
statistics in the antineutrino sample as opposed to a relatively high background. In 
this case the combination of neutrino and antineutrino fluxes is not as effective as for 
positive 5 and the allowed region increases (something that has been observed for the 
/3-Beam also, see Refs. [59,20]). 
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Fig. 16. Appearance and disappearance in the (9i‘i,6) plane at the Neutrino Factory. Left: 
L = 3000 km baseline; right: L = 7000 km baseline; top: 5 = 45°; bottom: 5 = —90°. The 
mixing angle takes the following values: @13 = 2°, 5° and 8 °. 


5 Sensitivities 


In this section we compare the sensitivities to ^ 13 , the maximality of 6 * 23 , the sign of 
the atmospheric mass difference, the 6 ' 23 -octant, as well as the CP-discovery potential 
potential at the experimental facilities studied in this paper. Since T2K-I and NOuA 
will run with neutrinos only, no sensitivity to the 6 phase is expected; we then study 
them separately from the other facilities in section (Sect.5.1). The SPL, T2K-II and the 
Neutrino Factories with L = 3000 km and L = 7000 km will be analysed in Sect.5.2. All 
the contours have been computed adding the appearance and disappearance channels 
available at a given facility and, if not differently stated, all the parameters not involved 
in the hts are hxed to their best £t values as quoted in the Introduction. 
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Fig. 17. Left: 3a 9insensitivity; right: 3a sensitivity to maximal 6 * 23 . Solid lines refer to 

T2K-I; dashed lines to NOuA. 

5.1 Sensitivities at T2K-I and NOuA 

In Fig. 17 we compare the 6 ' 13 -sensitivity (left plot) and the sensitivity to maximal 623 

(right plots) at T 2 K-I (solid lines) and NOz/A (dashed lines). 

• 6'13-sensitivity 

The 6'13-sensitivity is dehned as the one-parameter 3 a excluded region as a function 
of 6 in case of absence of signal. The contours presented in the following hgures 
represent the excluded valuees of 6'i3 for a given facility taking into account all 
posible choices of Satm and Soct- For both facilities the 6'13-sensitivity is rather poor, 
with a typical excluded region ranging from [sin^ Oislmin = [ 2 , 8] x 10 “^. Sensitivity 
is slightly better for negative 6 than for positive 6, with a maximal sensitivity at 
6 = — 90 °(— 135 °) for T 2 K (NOz/A). We have noticed that, while the T 2 K-I sensitivity 
is basically unaffected by the Satm, Soct choice, the NOz/A sensitivity is sensibly 
diminished in 5 G [— 90 °, 90 °] when choosing a wrong value for the sign of Am23 
or the 6'23-octant. 

To increase both the 6'13-sensitivity and the CP-discovery potential it has been 
proposed to add a proton driver to the NOz/A design, such as to increase the neu¬ 
trino flux from four to six times [ 26 , 40 ]. 

• Sensitivity to maximal 623 

The potential to exclude maximal 623 has been computed in the following way: 
for a given Am23, we look for the largest value of 623 for which the two-parameter 
3 a contours do not touch 623 = 45 °. We have considered both octants of 623 and we 
have found a behavior approximately symmetric of the sensitivity; we then display 
our results in the sin^ 623 variable up to the maximal mixing value of 0 . 5 . 

As we can see in Fig .17 (right), the sensitivity for the two experiments is essentially 
the same and strongly decreases for low values of Am^g. At the best £t point Am^g = 
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Fig. 18. Left: 3a Bi^-sensitivity; right: 3a CP-discovery potential. Solid lines refer to the 
L = 3000 km Neutrino Factory; dashed lines to the L = 7000 km Neutrino Factory; dotted 
lines to T2K-II; dot-dashed lines to the (standard) SPL. 

2.5 X 10“^ eV^ deviations as small as 14% of sin^6'23 from maximal mixing could be 
established at both facilities. Notice that these curves have been computed for a hxed 
013 = 0 (and Satm = +1) so that matter effects in disappearance probabilities and 
those from the appearance channel > z/g are completely negligible. We have also 
checked that for 0i3 close to the current bound, none of the shown results changes 
drastically and that everything is basically independent on the mass hierarchy. 


5.2 Sensitivities at T2K-II, SPL and Neutrino Factories 

• 013 -sensitivity and CP-discovery potential 

In Fig.18 (left) we present the 0i3-sensitivity (computed as explained in the pre¬ 
vious section) for T2K-II (dotted line), the SPL (dot-dashed line), the Neutrino 
Factory at L = 3000 (solid line) and L = 7000 (dashed line). As we can see, the NF 
at 3000 km shows the best sensitivity to 0i3 in the whole range of 5 (with the best 
value of sin^ 0i3 ~ 2 x 10“® for 5 around 30°) except a small region around 5 = —90° 
where T2K-II assures a slightly better performance (at the level of 8 x 10“®). For 
the Neutrino Factory at L = 7000 km, the loss in sensitivity is to be abscribed to 
the sign degeneracy. For example, we have found that for Satm = Satm we get that 0i3 
is excluded at 3a down to sin^ 0i3 ~ 7 x 10“® (0i3 ~ 0.5°) if 5 = 90°. On the other 
hand, for the same value of S (but for a wrong choice of Satm) we can only exclude 
013 down to sin^0i3 ~ 4 x 10““^ (0i3 ~ 1.1°). As for the SPL 0i3-sensitivity, the 3a 
excluded region varies in the range sin^ 0i3 G [4 x 10““^, 1.5 x 10“®] for different values 
of 6, with no big differences coming from different choices of Satm and Soct- T2K-II is 
signihcantly better than the SPL but signihcantly worse than the NF at 3000 km, 
apart from 6 = —90° (where it improves a little the NF limit). 

The CP-violation discovery potential has been computed as in Refs. [20,15]: at a 
hxed 013, we look for the smallest (largest) value of |5| for which the two-parameters 
3a contours of any of the degenerate solutions do not touch 5 = 0° nor |5| = 180°. No- 
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tice that, although the input 6^13 value is hxed, the degeneracies can touch 5 = 0 °, 180° 
at 9 i 3 7 ^ i 9 i 3 , also^. The outcome of this procedure is dually plotted, representing 
the region in the (^ 13 , 5) parameter space for which a CP-violating signal is observed 
at 3a. For each facility we show one single curve obtained taking into account the 
impact of all the degeneracies. 

First of all, notice that no CP-discovery potential has been considered for the 
NF at 7000 km: due to matter effects and the choice of the baseline (close to the 
magic baseline^ see Ref. [32]), the sensitivity to 5 vanishes. The best sensitivity to the 
measurement of the CP-violating phase is reached by the T2K-II experiment: a CP- 
violating signal can be observed at 3a for |(5| G [27°, 155°] down to sin^ 6^13 ~ 2 x 10“^ 
(^13 0.8°). At the SPL a good sensitivity can be reached also: CP violation is clearly 

observed at 3a for |5| E [45°, 135°] down to sin^ 6^13 ~ 9 x 10“^ ( 6 'i 3 ~ 1.7°), as it was 
shown in Ref. [20]. The NF at 3000 km is not as good as one would expect from its 
6 *i 3 -sensitivity contours, having a CP-discovery potential very similar to the SPL; the 
loss in sensitivity with respect to T2K-II is mainly due to the presence of the sign 
degeneracy. We have observed that the NF starts to be insensitive to the leptonic 
CP-violation just when, for ^13 ~ 2 °, sign degeneracies close to 6^13 = 0 appear, which 
do not allow to exclude the CP-conserving case any longer. This is a consequence 
of the fact that, for neutrinos with < >~ 30 GeV, the effective ^13 in matter 

[69,10] is extremely small. As it has been shown in Ref. [34], when 613 vanishes, the 
sign degeneracies flow to 5 = 0° or 5 = 180°, as can be seen in Figs. 2 and 3 (left) 
of Ref. [34]. At the SPL and T2K-II, for which the vacuum parameter 6*13 is the 
relevant parameter, this happens at a much smaller value ^13 ~ 0.5°. Therefore, at 
these experiments the CP-discovery potential is statistics-dominated. 

• Sensitivity to the mass hierarchy 

In addition to the previous sensitivities one can ask for the sensitivity to the sign 
of the atmospheric mass difference. We compute the smallest values of ^13 for which 
the sign of can be measured in the (sin^ 0 i 3 , 5) plane. The measurement of 

Satm requires matter effects to be sizeble. For this reason neither T2K-II nor the 
SPL have the capabilities to perform such a measurement. In Fig.19 we present the 
results for the Neutrino Factory only and for both posibilities of the true Satm = ±1. 
As expected, the NF at 3000 km (solid and dotted lines for normal and inverted 
hierarchy, respectively) exhibits the worst sensitivity to sign(Am 23 ). For Satm = +1 
the best sensitivity is reached in a quite large region around S = 90° for which 
[sin^ 6 *i 3 ]mm = 3.7 X 10“^ (1-1°) whereas for 6 ~ —90° [sin^ 6 *i 3 ]mm = 4 x 10“^ (3.6°). 
The situation is completely reversed for Satm = — 1 due to the fact that in matter a 


^ This is not the case of Fig. 11 in Ref. [59], where the excluded region in S at fixed 613 
in the absence of a CP-violating signal at 90% CL is presented. In practice, in that figure 
we compare N±{ 9 i 3 , with N±{ 6 i 3 , 0°), thus obtaining a one-parameter sensitivity plot in 5 
only. 
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Fig. 19. 3(7 sensitivity to the sign(Am 2 ^). Solid lines refer to the L = 3000 km Neutrino 
Factory with normal hierarchy; dashed lines to the L = 7000 km Neutrino Factory with 
normal hierarchy; dotted lines to the L = 3000 km Neutrino Factory with inverted hierarchy; 
dot-dashed lines to the L = 7000 km Neutrino Factory with inverted hierarchy. 

flip in the sign of Am^g corresponds to a change among neutrino and antineutrinos 
oscillation probabilities. At least one order of magnitude in sensitivity can be gained 
at the NF at 7000 km, depending on the fact that around L=7000 Km the effect of 
6 in the —>■ oscillation probability vanishes and that helps to measure the sign 
of the atmospheric mass difference. For the normal hierachy, a maximal sensitivity 
is achieved for 0 < 5 < 90° at the level of [sin^^igjmm = 4.4 x 10“^ (0.4°) whereas 
the largest value for [sin^ 0i3]min = 1-7 x 10“^ (0.8°) is around 6 = —100°. For the 
inverted hierarchy the sensitivity varies in the range sin^ 613 > [1, 3] x 10““^ (0.6°, 1°). 

• Sensitivity to maximal 623 and the octant-discovery potential 

Eventually in Fig.20 we present the sensitivity to maximal 623 (plot on the left) and 
the octant-discovery potential (plot on the right). The curves have been computed 
for 6*13 = 0 . 

With respect to the sensitivity to maximal 6*23, we observe that the NF at 3000 
km (solid line) is not as good as one may expect since it measures far away from 
the oscillation maximum. At the value of Am 23 in which the sensitivity is maximal, 
deviations as small as 10 % of sin^ 6*23 from maximal mixing could be established. 
Similar behaviour is expected at the SPL (dot-dashed line), except for small Am^g in 
which it outperforms the NF at 3000 km. A big improvement will be however achieved 
at a NF at 7000 km (dashed lines) and T2K-II (dotted line); both experiments have 
energies and baselines (as well as off-axis angle for T2K-II) chosen to match the hrst 
oscillation peak (in vacuum and matter respectively); the reached sensitivities are at 
the level of sin^ 6*23 G [0.45 — 0.48] almost independently on the value of Am^g, which 
means that deviations from maximal mixing of the order of 4% could be established. 

Notice that, although this sensitivity is rather good, in general it is very difficult to 
determine the octant in which the atmospheric angle lies. As we can see from eq.(l), 
it is quite difficult to break the 6*23 —>■ n/2 — 623 symmetry induced by the leading 
term in the transition probability; the subleading terms that could help in lifting this 
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Fig. 20. Left: 3a sensitivity to maximal 023/ right: 3a sensitivity to the 823 -octant. Solid lines 
refer to the L = 3000 km Neutrino Factory; dashed lines to the L = 7000 km Neutrino 
Factory; dotted lines to T2K-II; dot-dashed lines to the (standard) SPL. 

degeneracy are very difficult to isolate. However, for values of 613 different from zero, 
we can take full advantage of matter effects in the disappearance of muon neutrinos, 
as we have already seen in Sect.4 (Figs. 12-13). Obviously the matter effects at SPL 
and T2K can never be sizeble enough to solve the octant ambiguity; on the other 
hand, the Neutrino Factory shows a (limited) capability to solve it, irrespective of 
the baseline and the value of 6. To illustrate this point, we hxed 623 = 41.5° and 
= 2.5 x 10“^ eV^ and, for any value of S, we compute the minimum value 
of 013 for which the octant ambiguity is solved. As we can see in the right plot of 
Fig.20, a 3a octant discovery is possible for sin^ 023 > 0.01 (6°). 


6 Conclusions 

In this paper, we studied the measurement of the atmospheric neutrino oscillation 
parameters, 623 and Am^a using the disappearance channel at three conventional 
beam facilities, the SPL, T2K-phase I and NOi/A, and at the Neutrino Factory. The 
precision on these two parameters will be of crucial importance in the measurement of 
two of the unknowns of the PMNS mixing matrix, 0i3 and the leptonic CP-violating 
phase S. 

It has been shown that counting experiments cannot reduce significantly the uncer¬ 
tainties on 023 and Am^a. Hence, we have considered detectors with (modest but non 
vanishing) sensitivity to the energy of the final leptons produced via neutrino interac¬ 
tion. For T2K-I and NOz/A we have indeed found that, independently of the input value 
of 023, the errors on 623 can be signihcantly reduced. Different energy bins give different 
allowed regions in the ( 023 , Am 23 )-plane and the combination of them eventually helps 
in reducing the uncertainties around the physical point. The two experiments only run 
with TT’*', and thus are not sensitive to the CP-violating phase 6 . The SPL, with a 1 
Mton detector and 2-1-8 years of running time with tt’*' and 7 i~, respectively, could 
in principle greatly improve these results. However, having neutrinos and antineutri- 
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nos at this setup does not help, since they have been produced with the same energy, 
roughly. Moreover, the baseline is too short for matter effects to make a difference. 
As a consequence, information coming from the two fluxes is not complementary and 
just add statistically. Remember, however, that this is not the case in the appearance 
mode, z/g, for which this setup was designed. In this paper, we have considered 

two bins of energy at the SPL: this already makes a big difference in the precision on 
023 as it can be seen comparing with our results of Refs. [59,20] and [15], where only a 
counting experiment was considered for this facility. 

We have then studied the impact of ^13 and S in the measurement of 6*23 and Aml^. 
Since both parameters are unknown at present, their values should be reconstructed at 
the same time with the atmospheric parameters. We have shown that the main effect 
of considering ^13 as a free variable is a shift of the contours toward larger values of 
023 (with respect to the input point) while a free 6 in the whole 6 G [—180°, 180°] 
range does not produce any significant distortion of the allowed regions. These effects 
can be strongly softened if, in addition to the disappearance channels, we introduce 
in the analysis the appearance z/^ —z/e channels, which the Super-Beam facilities 
have been proposed to look for. Two important conclusions can then be drawn: if 6^13 
is kept fixed, the combination of appearance and disappearance channels solves the 
disappearance octant degeneracy for 6*13 large enough; on the other hand, the octant 
degeneracy is not lifted if 0 i 3 is free to vary in the current allowed range, for almost 
any value of 6 ^ 13 . 

The situation is quite different at the Neutrino Factory. We have pointed out that the 
disappearance channel is a very powerful tool to reduce the uncertainties on 
the atmospheric parameters up to an unprecedented level and, in combination with 
appearance channels, to solve many of the discrete ambiguities affecting the measure¬ 
ment of the PMNS matrix elements. The main feature emerging from such an analysis 
is that, for ^13 > 3° — 4°, the sinergy between disappearance and appearance channels 
(the more renowned golden and silver channels) greatly helps in solving the octant 
ambiguity, at both L = 3000 km and L = 7000 km baselines, as a result of the strong 
matter effects along the neutrino path and the huge statistics at hand. This remains 
true even in the case we leave ^13 as a free parameter, a situation which has been shown 
not to be true for the other facilities. Equally remarkable, in the same range of ^ 13 , the 
sign clones are solved and they are not present in the fits, independently of the baseline 
and the input values of 023 and This allows a precise measurement of 0 i3 at both 

facilities while the CP-phase S, as it is already well known, can only be measured at 
L = 3000 Km. For values of 0 i3 < 3°, some remnant of the clones still remains and the 
measurement of the two unknowns cannot be performed with huge precision. 

Eventually, we have presented a comparison at the different facilities of the 3a sen¬ 
sitivity to 013 , to maximal 023 , to the mass hierarchy, to the 6 ' 23 -octant and their 
CP-discovery potential. T2K-I and NOz/A, running with neutrinos only are not able 
to measure 6 . On the other hand, they have very similar performances in the 013- 
sensitivity, with an excluded region ranging from [sin^ 013]min = 2 X 10 3 to 8 X 10 ^ 
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and in the sensitivity to maximal 6 ^ 23 , with the maximal ^23 that can be distinguished 
from 6*23 = 45° ranging from [sin^ 6 * 23 ]max = 0.45 to 0.40 for G [2.0, 3.0] x 10“^ 

eV2. 

The Neutrino Factory outperforms the considered Super-Beams for both baselines in 
the 6 *i 3 -sensitivity and in the sensitivity to the sign of the atmospheric mass difference. 
The longest baseline is as effective as T2K-II in the sensitivity to maximal 6 * 23 . The 
excluded region in 6*13 at the 3000 km Neutrino Factory ranges from [sin^ 6 * 13 ]mm = 
2 X 10“^ to 2 X 10“^ (for particular values of 5). The sign of the atmospheric mass 
difference can be measured at the 7000 Km (3000 Km) Neutrino Factory for sin^ 6*13 
as small as 10“"^ (10“^). Finally, the 023-octant can be identihed at both facilities for 
sin^ 0 i 3 as small as 10“^ ( 0 i 3 > 6 °). Notably enough, the CP-discovery potential seems 
to be larger for T2K-II than for Neutrino Factory at L = 3000 Km (no CP-discovery 
potential is expected at the magic baseline L = 7000 Km). This is explained as follows: 
as a general rule, for small values of 0 i 3 the degeneracies flow toward 5 = 0° and 
|5| = 180° (see Ref. [34] and Ref. [33]), thus mimicking a non-CP violating phase. Due 
to a parametric conspiracy between the chosen energy and baseline and the matter 
effects, at the Neutrino Factory the typical value of 0 i 3 for which this happens is much 
larger than at the SPL and T2K. Therefore, although from the statistical point of 
view the Neutrino Factory would certainly outperform both the SPL and T2K-II, in 
practice for small values of 0 i 3 a CP-violating phase will be difficult to distinguish from 
a non-CP-violating one, if Satm and Soct are not measured previously. 


Acknowledgments 

We would like to thank E. Couce, B. Gavela, J. Gomez-Cadenas, P. Hernandez, P. Hu¬ 
ber, O. Mena, P. Migliozzi, T. Schwetz and W. Winter for useful discussions. The 
authors acknowledge the financial support of MCYT through project FPA2003-04597 
and of the European Union through the networking activity BENE. E.F. acknowledges 
hnancial support from the UAM. 


Appendix A 

In this Appendix we recall the number of events per bin that are expected in the 
disappearance channel at the SPL, the T2K-I, the NOz/A and the Neutrino 
Factory experiments. The expected T2K-I muon and electron identihcation efficiencies 
are presented, as well. 

In all tables, two different values of 0 i 3 and S are considered, as well as two choices of 
sign(Am|3). 

In Tab. 1 we show the expected event rates per bin at the SPL after a 2 years run with 
and a 8 years run with 7 r“, using the fluxes of Ref. [43] and a 440 kt detector at 
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L = 130 km. 
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Table 1 

Disappearance event rates for a 2 + 8 years vr"'' + tt~ run at the standard SPL, for different 
values of 013,5 and of the sign of the atmospheric mass difference, Satm- 

In Tab. 2 we show the expected event rates for the two lowest bins at T2K-I after a 5 
years run with 7 r~^, for a 22.5 kt detector at L = 295 km. 
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Table 2 

Disappearance event rates for a 5 years 7r+ run at T2K-I, for different values 0/6*13, ^ o,nd of 
the sign of the atmospheric mass difference, Satm- 

In Tab. 3 we show the expected event rates per bin at NOz/A after a 5 years run with 
Ti^, for a 30 kt detector at L = 812 km. 

In Tabs. 4 and 5 we show the expected event rates per bin after a 5 years run with 
and a 5 years run with pff at the Neutrino Factory, for a 40 kt detector at L = 3000 
and L = 7000 km, respectively. 
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Table 3 

Disappearance event rates for a 5 years run at NOvA, for different values of 613,6 and of 
the sign of the atmospheric mass difference, Satm- 
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Table 4 

Disappearance event rates for a 5 + 5 years run at the L = 3000 km Neutrino Factory, for 
different values of 613,6 and of the sign of the atmospheric mass difference, Satm- 


Eventually, in Fig. 21 we present the muon (solid) and electron (dashed) identihcation 
efficiencies at the T2K-I experiment, as found in Ref. [70] and [16], respectively. 
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Table 5 

Disappearance event rates for a 5 + 5 years run at the L = 7000 km Neutrino Factory, for 
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Fig. 21. Muon (solid) and electron (dashed) identification efficiency at T2K-I as a function 
of the neutrino energy. 
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